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Abstract 

We consider Galileon models on curved spacetime, as well as the counterterms introduced to 
maintain the second-order nature of the field equations of these models when both the metric and 
the scalar are made dynamical. Working in a gauge invariant framework, we first show how all 
the third-order time derivatives appearing in the field equations — both metric and scalar — of 
a Galileon model or one defined by a given counterterm can be eliminated to leave held equations 
which contain at most second-order time derivatives of the metric and of the scalar. The same 
is shown to hold for arbitrary linear combinations of such models, as well as their k-essence- 
like/Horndeski generalizations. This supports the claim that the number of degrees of freedom 
in these models is only 3, counting 2 for the graviton and 1 for the scalar. We comment on the 
arguments given previously in support of this claim. We then prove that this number of degrees of 
freedom is strictly less that 4 in one particular such model by carrying out a full-hedged Hamiltonian 
analysis. In contrast to previous results, our analyses do not assume any particular gauge choice 
of restricted applicability. 
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I. INTRODUCTION 


Scalar-tensor theories are widely used in cosmology and extensions of general relativity, 
with applications ranging from inflation to the late-time observed acceleration of the Uni¬ 
verse, and tests of gravitation. Motivated in part by the ability of some of these theories 
to give an alternative to dark energy, there has recently been renewed interest in the delin¬ 
eation of theories in which — besides the two degrees of freedom of a massless graviton — 
there is only one propagating degree of freedom (d.o.f.) stemming from the scalar. Along 
this line, an important result was achieved by Horndeski who classified all scalar-tensor 
theories in 4 dimensions having field equations (both for the metric and the scalar) with 
derivatives of order less than or equal to two [I|. Similarly, Ref. (see also for earlier 
works) introduced, on flat spacetime and for an arbitrary number of dimensions D, a set 
of scalar theories with field equations exactly of second order: the Galileons. These the¬ 
ories were later “covariantized”, i.e., put on arbitrary curved spacetime with a dynamical 
metric, while maintaining the second-order nature of the scalar field equation, as well as 
yielding metric field equations of the same order 1,0. This covariantization procedure 
involves a non-minimal coupling between the curvature and the scalar in the form of very 
specihc counterterms able to remove all higher derivatives from the held equations. Indeed 
a minimal covariantization of the original Galileon of Ref. 0 (i.e., the mere replacement of 
partial derivatives by covariant derivatives in the action) was shown to lead to third-order 
derivatives in the held equations.^ Another relevant work is that of 0, which classihed 
all scalar theories having equations of motion of order less than or equal to two on a hat 
spacetime of arbitrary dimension, and then covariantized these theories. It was shown that 
the original hat-space time Galileons 0, their hat spacetime generalizations 0 as well as 
their covariantization 00 (with the meaning above) belong, for a spacetime with 4 dimen¬ 
sions, to the set of Horndeski (as they should according to Horndeski theorem) 0,0. These 
theories were also generalized to the case of multiscalars and p-forms @-[0 


Having covariant second-order held equations is a priori enough, once diheomorphism 
invariance is taken into account, to have just 3 propagating degrees of freedom in vacuum 



minimally covariantized Galileons, the number of propagating degrees of freedom can still 
be only three due to the presence of some hidden constraint in the theory. (In minimally 
covariantized Galileons and related models, as stressed in 0, 0 , the held equations for the 


^ We refer to such models as “minimally covariantized” Galileons, to stress the difference with the “covari¬ 
antized” Galileons of ref. [l|, 0, which contain non-minimal coupling to curvature in the form of the 
counterterms mentioned above. 
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metric are second order for the metric but contain third-order derivatives of the scalar, and 
conversely the scalar held equation is second order for the scalar but third-order for the 
metric.) This appears, of course, to be a perfectly legitimate possibility and it is not hard 
to build some simple examples with such a property (see e.g. 12811) . However, so far, the 
arguments given in favor of this happening, as proposed in Refs. |2l|, [22|, \2j \ , do not appear 
to us to be entirely conclusive. 

The reasons are the following. First, the Hamiltonian analysis of the corresponding theory 


has only been carried starting from a gauge hxed action [2lH23| where (i) the scalar gradient 
Vfj,ip is assumed to give the time how direction, and (ii) this gauge hxing is not included 
in the Hamiltonian (i.e., the gauge is explicitly broken to start with). This gauge, usually 
referred to as the “unitary gauge”, hides all the dynamics of the scalar and it is easy to see 
that it eliminates all third time derivatives in the held equations (see Refs. 0,0 and also 
Sec. mil belowh In this sense, it is perhaps not surprising that working in this gauge, one 
hnds less degrees of freedom than those expected from an Ostrogradski-type of reasoning. 
Furthermore, this unitary gauge choice is obviously only possible if the scalar gradient is 
everywhere time-like (or at least time-like in the vicinity of some would-be Cauchy surface), 
a situation which only covers rather limited subset of all possibilities. Indeed, it does not 
allow one to say anything about the Cauchy problem when, on the Cauchy surface, the 
scalar has a gradient which is not always time-like — a perfectly legitimate choice of initial 
condition. For instance, a physical situation of major importance of this kind is that of a 
static and spherically symmetric background, since the hypersurfaces (p = const, are not 
spacelike and cannot be chosen as initial value surfaces. Second, it is well known, e.g. when 
considering Maxwell theory, that analyzing the d.o.f. content of the theory in a given gauge 
can be very misleading, in particular when the gauge is explicitly broken to start with. 
Finally, the covariant reasoning given in (analogous to the one of jl^) appears to us to 
be incomplete if not incorrect. Indeed, there it is stated that taking an appropriate trace 
of the metric held equations (which are known to contain third-order time derivatives of 
the scalar) enables one to extract the third time derivatives of the scalar in terms of second 
derivatives — and then, inserting this back in the metric held equations, gives a second- 
order system (and similarly with the scalar held equation). This seems to omit the fact 
that in this way one can at best eliminate from the held equations all but one third-order 
time derivatives of the scalar (and similarly for the metric): the reason is that the trace 
of the metric held equations is itself a held equation which must still be solved, and which 
still contains a third-order time derivative. Hence, in contrast to the claims in j^. 27|. 
the covariant procedure outlined in those papers appears not to lead to a complete set of 
equations in which all third-order time derivatives have been eliminated. 

Here we will reexamine these issues and argue, in two diherent ways, that minimally 
covariantized Galileons indeed propagate less degrees of freedom than expected from the 
third-order nature of the held equations. Throughout we work in a totally gauge-invariant 
framework. This paper is organized as follows. In Sec. Hi] we show how the system of 
-|- 1 field equations^ of the theory considered — namely all minimally covariantized 
Galileons and independently all the counterterms, as well as any linear combination of them 
and their Horndeski-like generalizations — can be reduced to an equivalent system with only 
second-order time derivatives (using, however, a very diherent procedure from the one given 


^ This number -|-1 will be quoted several times in the present paper. In D = A dimensions, it simply 

reduces to the usual 10 -I- 1 = 11 field equations of the metric and the scalar field. 
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in Refs. j^ . 2^). In Sec. IIIIl we present a Hamiltonian analysis of one particular theory 
in the minimally covariantized Galileon family, namely the simplest non trivial one (in the 
sense that its held equations do contain third-order derivatives), to show that the number 
of constraints is sufficient to yield less than 4 propagating degrees of freedom. A last section 
gives our conclusions. 


II. REMOVAL OF THIRD TIME DERIVATIVES IN THE GENERAL CASE 


In a spacetime with metric in any dimension D, Ref. introduced the generalized 
Galileon Lagrangians C(n+i,p)- These involve a total of 2n derivatives acting on a product 
of n -|- 1 scalar helds (p, and p Riemann tensors R\^up- The Lagrangians with p = 0, the 
are the “minimally covariantized” Galileons, whereas the Lagrangians with p ^ 0 
are called the “counterterms”. 

For a given n, it was shown that — up to an irrelevant global factor — there exists a 
unique linear combination 


Pmax 

E 

p=0 


a 


{n+l,p)R{n+l,p) 


such that all held equations are of second order. These are the 
Here 

n — 1 

Pmax 


( 1 ) 


‘covariantized” Galileons. 


( 2 ) 


is the integer part of and the constant coefficients C(^n+i,p) take a very specihc form 
which may be found in Eq. (37) of [^. 

Gonversely, any other linear combination (for instance, each of these Lagrangians C(^n+i,p) 
individually when Pmax > 0) does yield third derivatives in its held equations. More specih- 
cally, the scalar held equation contains third time derivatives of the metric tensor, and the 
Einstein equations contain the third time derivative of the scalar held. To be able to com¬ 
pute the time evolution, it seems thus necessary to specify more initial data on a Gauchy 
surface, and one expects the existence of more degrees of freedom than just a single scalar 
held and the two helicities of the graviton. 

Although higher-order held equations indeed generically lead to extra degrees of freedom 


(which are even generically ghost modes jlTI. Il8| . implying the instability of the theory), 
specihc examples show that this is not always the case. There may, for instance, exist 
extra constraints (related or not to some hidden gauge symmetry) which kill some of the 
modes. Or the few equations involving third (or higher) time derivatives may actually be 
obtained by diherentiating, with respect to time, some independently known second-order 
held equations. In this case, extra initial data are not necessarily needed on the Gauchy 
surface. An elegant toy-model of this kind was presented in Sec. 7.1 of 
result applies to “mimetic dark matter’ 


28 , and a similar 


29 


The aim of this section is to show that the held equations of the generalized Galileon 
Lagrangians C(^n+i,p) are of this second kind: all the third time derivatives can be obtained 
by deriving independently known second-order held equations, and can therefore be re¬ 
moved from all held equations. This supports the main claims of Refs. j21-2^|, although 
our procedure dihers from theirs, and does not suher from the problems mentioned in the 
Introduction. In particular, we will not hx any gauge in our derivation. 
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Thus we consider the D-dimensional theories defined by the action 

S = (Seh + ^Gal (3) 

in curved spacetime, where Seh = / d^x^J—gR is the Einstein-Hilbert action^, without any 
factor c^/lbvrG to simplify our discussion, and where 


5'Gal = 


d^xy/—g 



(4) 


with arbitrary constant coefficients fc(n,p) — and hence not the specific C(n,p) discussed above. 
(We will furthermore consider the Horndeski-like generalization of these theories at the end 
of the section.) We also define the Galileon energy-momentum tensor as 


_ 


1 5S, 


Gal 




flU 


(5) 


without any factor 2, so that Einstein’s equations (i.e., the field equations for the metric) 
simply read where is the Einstein tensor. Finally, we define £ = SScai/dip, 

so that the scalar field equation reads £ = 0. 

Note that this scalar field equation is a consequence of Einstein’s equations since, because 
of the diffeomorphism invariance of action S, it follows that 


-T'^^), (6) 

where (p'^ = denotes the covariant derivative of the scalar field (without writing any 
semicolon, to simplify; we shall also write = V... in the following). Indepen¬ 

dently of the diffeomorphism-invariance argument, this can also be checked explicitly for 
the general Lagrangians C{n,p) or for particular cases [^. Therefore, if we manage to prove 
that Einstein’s equations can be recast as a set of second-order differential equations (with 
respect to time), Eq. ([6]) shows that the third time derivatives of the metric tensor entering 
the scalar field equation £ = 0 should not pose more problems than in the toy-model of 
Ref. [ 2 ^. It should be noted that at the spacetime points where p'^ happens to vanish, £ 
can no longer be extracted from (jS]). However, it is easy to see that in all field equations, 
third derivatives are always multiplied by a gradient p'' (and even several of them). There¬ 
fore, at the points where p^ = 0, all field equations are at most of second order, and do not 
pose any problem. In addition to the above argument based on Eq. ([6]), we will actually 
prove an even stronger result below: It is also possible to recast the scalar field equation 
itself as a second-order one (with respect to time), by combining it with the time derivative 
of another linear combination of Einstein’s equations. 


A. Example: in If = 4 dimensions 

Before attacking the general case of £(n,p) which is rather technically involved, we begin 
in this subsection by illustrating how the steps work in D = 4 dimensions, focusing on the 
simplest non-trivial Galileon action, namely S'Gai = / d'^xy^—gC(4^oy 


^ Throughout this paper, we use the sign conventions of Ref. 


30|, and notably the mostly-plus signature. 
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The Lagrangian >C( 4 ,o); which is also the subject of the Hamiltonian analysis of Sec. IIIIl 
is given by 


£(4,0) = - 2 (□</9) (v?pV?'’) + 2 




(7) 


where is the Levi-Civita (fully antisymmetric) tensor in 4 dimensions [see Eq. fl2S]l for 
its definition in D dimensions]. Its stress-energy tensor is given by 


rjiflU 

^ (4,0) 


pi'^piP3P5a^/i2A‘4M6 


-q^ e 










-£pi 




( 8 ) 


In the remainder of this subsection we simply denote £( 4 ’))) by 

The first step is to determine explicitly the different order time derivatives which appear 
in the zero-zero component of the Einstein equation recalling that in there 

are at most first time derivatives of the metric. From ([ 8 ]) it follows that 


from which we immediately see that contains no terms in V? nor in (pi (where latin 
indices mean spatial components). This latter term could be generated from the term within 
square brackets, but that would require /is = p-s = 0 , in which case the result vanishes by 
antisymmetry of the Levi-Civita tensor. Thus the highest order time derivative of the scalar 
field it contains is (/3, whose coefficient can be determined directly from (|9]) and will be 
given below. Regarding the metric, there are second-order time derivatives coming from 
the third-order covariant derivatives of tp on the second line, since b) We 

must therefore take pa = 0 and a = 0 to find these, and we obtain 

= w . ( 10 ) 


where N = is the usual lapse in the ADM decomposition, see Eq. flAip . and 

is the 3-dimensional Levi-Civita tensor related to the 4-dimensional one by 

ijk 

( 11 ) 

The subscript gij on the left hand side of (1T0|) is due to the fact that the Christoffel symbols 
Tgj- only contain first time derivatives of the spatial components of the metric (see Appendix 
More explicitly 

{doT'^j)pi, = —Np^{doKqj) + first-order derivatives, ( 12 ) 


^ Hence note that contains no terms in N or Ni (where N and M are the usual lapse and shift in the 
ADM decomposition). Hence there are no terms in ijoo nor goi- 
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where Kij is the extrinsic curvature. Thus 

^"ls« = . (13) 

The second step involves carrying out the same procedure for (piT^\ From ([8]) with p = 0 
and u = i we hud 

-'Pn'Pi - •Pn'Pi [£'"'‘’'‘'‘£'“''t('9'Wrea);w] 

+ 9 ,,^"’’'=^'“'“ tiWW,, ‘P,„ ):,■ ( 14 ) 

Again, it is clear that there are no terms in ‘F. Similarly there are no terms in ipj since one 
always ends up with a contraction = 0. There are obviously terms in ip (see below). 

Concerning the terms in Qij, following the same logic as above, we hnd that they are given 
by 

^ . (16) 

It follows from flTOll and flTHll that and contain exactly the same combination of 

second-order derivatives of the metric. 

Furthermore, we recall that the components of the Einstein tensor do not involve 
second time derivatives of the metric. Indeed, the Bianchi identities imply the covariant 
conservation of the Einstein tensor, = 0, therefore 

^0^°^ = -diG^>^ - 0(FG), (16) 

where 0{TG) means the four terms involving contractions of Christoffel symbols F^^ with 
the Einstein tensor G^^. Since the right-hand side of Eq. ffT6l) contains at most second time 
derivatives of the metric tensor, this must be so for the left-hand side doG^^, therefore the 
components G^^ contain at most hrst time derivatives. 

Thus we arrive at the hrst important conclusion that, using Eqs. oa and fllSp ' th.6 
combination of the Einstein equations 

- ((p*(p‘)T°° = = p\Pt.G^n - {P.P^G^" (17) 

determines p in terms of first time derivatives. As a result, any P appearing in the equations 
of motion can be expressed in terms of second time derivatives simply by differentiating flTT)) . 

The next step is the determination of the coefficients of these p terms. Starting from ([9]) 
and flTT)) we hnd 

^ kPiPmidoT'^j)pfi\ + Bpo{p^p) + lower order derivatives, (18) 

p^T^t^ = [£"''"^£'hPiPm{doT'^j)pu] - A{P^P) + lower order derivatives, (19) 

where 

® . 61 = ( 20 ) 


7 






(Note that A is in fact nothing other than >C( 3 ,o) with only spatial indices.) Hence the second 
combination of Einstein’s equations 

( 21 ) 

is an equation for the specihc combination of gij appearing in 

( 22 ) 

in terms of hrst time derivatives of the helds. 

As a hnal step we must show that the third-order time derivatives of the metric appearing 
in the equation of motion £ = Q for the scalar held, are exactly given by the time derivative 
of the combination appearing in ([22]). The third-order derivatives of the metric in £ are 

s ~ (23) 

where Rij,upa is the Riemann tensor. To hnd the 9 ij terms appearing here, it is sufficient to 
set A = 0, and then the relevant third-order time derivative is simply the derivative of the 
term in g appearing in 

• (24) 

However, using the results of Appendix [^ this is nothing other than the combination fl22|) 
(up to irrelevant numerical factors). Since this is the contraction appearing in fl2Tl) . it can 
thus be expressed in terms of hrst time derivatives of the held. Thus the third-order time 
derivatives of the metric appearing in the scalar held equation of motion can be replaced by 
second-order time derivatives on using the derivative of fl^TD . 

Hence we arrive at the conclusion that, despite containing higher-order time derivatives, 
all 11 equations of motion for this theory can be expressed solely in terms of second-order 
time derivatives: the derivative of the combination ffT7|) gives V? in terms of second-order 
time derivatives, whilst the derivative of the combination (|2T|) gives the required contraction 
of g ij in terms of second-order time derivatives. 

We now generalize these results to arbitrary D, n and p, which requires us to intro¬ 
duce more powerful notation. We will also consider arbitrary linear combinations of these 
Lagrangians. 


B. General case: arbitrary D, n and p 


For the following discussion, the most convenient expression for the Lagrangians 
C(n+i,p) uses the Levi-Civita fully antisymmetric tensor 




1 

\/^ 


[Mi C-/42 




..6 


Md] 
D ) 


(25) 


where the square bracket denotes unnormalized permutations. In any dimension D > n, we 
dehne 

B{n+l,p) £^{2n)^fj,i^p2'^ip)^{q)^ (26) 

where A( 2 n) is a compact notation for 


•^(2n) 


_A‘lM3A‘5---M2n-l l^ll'2-"^D-n M2M4M6■ ■ ■M277 


{D-n)\ 


( 27 ) 






and where 


"^(p) 

‘ 5 ( 9 ) 




V-ii-1 P4i + 1 P4i M4i+2) 


i=l 


9-1 

n 2i ^^2n-2i 1 

i=0 


with 

g = n — 1 — 2 p. 

These dehnitions assume that n, p and q are integers in the ranges 

n — 1 


1 < n < D, 1 < p < 


, and 1 < g < n — 1. 


(28) 

(29) 


(30) 


(31) 


We also set 7^(o) = 1 for p = 0 and iS(o) = 1 for g = 0, and use the convention that 7?.(p) = 0 
for p < 0 and iS(g) = 0 for g < 0. [The cases n = 0 and n = —1, with p = 0, may also 
be dehned as = (p and /l(o,o) = const., but we do not consider them here since 

they obviously do not yield higher-order held equations.] The numerical factor 1/{D — n)\ 
is introduced in Eq. fl27D so that A{ 2 n) keeps the same expression in terms of products of 
metric tensors in any dimension. 

When p = 0, i.e., without any Riemann tensor involved, these dehnitions reduce to the 
Galileons of Ref. [^. For instance, /l( 2 ,o) = is the kinetic term of a standard scalar held 
(though, referring to Eq. (jl]), one should choose a negative k{ 2 fl) in order not to have a ghost 



valid for all C[n+i,p) in any dimension D, as well as linear combinations of them. 


1 . C[n+i,p) mid their linear combinations 

We hrst focus on a single Lagrangian C(^n+i,p)- In order to write its Einstein equations 
Qa/3 _ j^a/3 simplest way, it will be useful to introduce the following compact notation 

generalizing fl27|) : 




_ Anin2---ni-ini+i---n2n 


= A‘ 

— •^(2n) ’ 


(2n) 


(32) 

(33) 
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where 2 n is the rank of the tensor, and i and j locate the positions of the indices a and f5 
which are explicitly indicated. The energy-momentum tensor is then given by 


rjnOi0 _ 


2 V?! 


^{2n) — 




(2n+2,2n+l,2n+2) 

'TX / 

■*"2'^{2^-4p+3,4p+4) \}P V^Pl‘/^P 2 ”^(p)‘^{g-l)] ;a- 

-|-2P'^(2n,4p+l,4p+2) . 

—pA 


^Pi^P 2 "^(p)‘^{g) 


/J-dp f-^Ap—1 


(" uP) 

(2n,4p—1) P4p+1 P4p p4p+2 


9^P1^P2 (9^a)"^(p-1)‘^(ij)) 


(34) 


where symmetrization over a and is assumed [i.e., = (X“^-|-X ^“)/2 for any tensor], 

notably on the second and last two lines which are not automatically symmetric. Regarding 
the first term, note that the factor p/p\ does not cause any divergence: it vanishes when 
p = 0, and when p 7 ^ 0 it is multiplied by (<p|)^ contained in the Tt{p) term fl28|l . The tensor 

'^^2n+2 2n+i 2n+2) 2 ? 7 , -|- 2 free iudices, the last two of them being a and 

(3. The corresponding term in originates from the variation with respect to the metric 
of the contracted indices in Eq. (HTh . It vanishes if D < n -|- 1. Note that in the second line, 
the factor has been extracted from the square bracket which is covariantly derived with 
respect to p^p+s, because ^ 4 ^ 4,3 with two odd indices would vanish when contracted with 
the hrst antisymmetric e tensor of Eq. fl27ll . 

This energy-momentum tensor in Eq. fIMl) allows us to prove several important lemmas. 
Let us focus to start with on third derivatives: by suitably permuting and relabeling dummy 
indices, one may rewrite them as 


rpa(3 

3rd der. 


9(9 - 1 ) J^aP ^ c ,.A,„ 

2 ‘^(2n,4p-|-3,4p+4)^Pl^P2 'Mp)‘^(<?- 2)9^ ‘^p4p+5 P4p+6 A 

-|-4p •A(2n,4p+l,4p-|-2)^Ml^P2'^(p-l)‘^(g)V^ V^Ap4pP4p-l- 


(35) 


This is is equivalent to Eq. (34) of Ref. . We can thus arrive at the following conclusions. 

First, no third-order derivative of the metric tensor enters T"^. Indeed, the Bianchi 
identities R\p,[yp-,a] = 0 cancel most of the differentiated Riemann tensors in fl341) . because 
three of their indices are contracted with the same antisymmetric e tensor of . 4 ,( 2 n)- The 
only non vanishing ones come from the 'R-{p)-^a of the third term, but they exactly cancel with 
the derivatives of the Riemann tensors generated by permuting the indices of 
coming from the fourth term. 

Second, Eq. fl35|) also proves that does not contain any <P nor p*. Indeed, because 
CK = /3 = 0, one odd index (4p -|- 3 or 4p -|- 1) and one even index (4p -|- 4 or 4p -|- 2) must 
be 0 , therefore all other (contracted) indices of the two e tensors in A{ 2 n) must be spatial. 
In conclusion, the only third-differentiated scalar fields fl35l) cannot contain more than one 
time derivative in 

Third, the contraction (and in particular (po,T“°, that we will use in the argument 

below) does not contain any third derivative of the scalar held. Indeed, ipa and are 
then contracted with the same antisymmetric £ tensor entering A( 2 n), m Eq. fl35p . making 
it vanish. 

We now proceed as in Section III Al and compute the second time derivatives entering 
<PaT"° and T°°, using the full expression fl3T|) . However, before doing so, and in order to 
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simplify the resulting expressions, we introduce the following further notation. Using ([33]), 
we dehne 

'^(n+l,p) = ~•^^n+2,2n+l,2n+2)‘^^ll‘^^^2^{p)^{,q)1 ( 36 ) 

which is the analogue of Eq. fl26p . but where now all contracted indices are spatial (though 
the covariant derivatives and the Riemann tensors remain D-dimensional). As before, see 
Eq. fl30|l . q = n — 1 — 2p. Note that A^2n+2 2n+i 2n+2) extra indices relative to A[2n)-i 

namely the last two which are both 0, meaning that the 2n Erst must be spatial. Similarly, 
we define 

^(n+l,p) — ~'^(2n+4,2n+3,2n+4) M 2 n +2 ’ (3^) 

which, as before, contains n-|-l scalar helds and p Riemann tensors, as well as now a Christof- 
fel symbol r°-. Notice that all contracted indices again are spatial, because the (2n -|- 3)rd 
and (2n-|-4)th indices are imposed to be 0. [Actually, contains —(n — ^y] 

On using flMD . the second time derivatives entering and can then be written 

as 


where 


V?a^2"nd t der. = “V?” {A(p + plC), 

rZtder. = 


A 

B 

C 


3) ^spatial I ^2/’Spatial 


(n,p—1) ’ 


?(?-l)/’r 

-L.(. 


2 ~'(n—l,p) A 4p 

^00 c ;; 

^ ■^(2n,4p+3,4p+4) Apt '^(p)‘^(g-2)5'p4p+5 M4p+6 

A2p •A.(-2/j^4p_1^4p)9^pi9^p2'^(p—l)‘^(g)(?M4p+l /44p+2 I 


(38) 

(39) 

(40) 

(41) 

(42) 


with g = n — 1 — 2p in these expressions — but the g’s involved within and depend 

on their precise indices (n,p) as dehned in Eqs. (136|) and (l37|) . Note that A and B involve at 
most hrst time derivatives, whereas C contains a specihc contraction of gij with other helds 
(themselves differentiated at most once with respect to time). 

The fact that the same contraction C of gij enters both fl38ll and fl39ll will be crucial below. 
Whilst this looks miraculous, it is actually simply a consequence of the diffeomorphism 
invariance of the theory, together with the fact that T“^ involves at most second derivatives 
of the metric as proved in Eq. fl35|) . Indeed, diffeomorphism invariance implies Eq. (|6]), 
which tells us that (p°'8 = 2doT^°' + 2cIjT*“ -|- 0{TT). The third time derivatives of the 
metric entering £ can thus come only from 2 doT^^, proving that all must contain 
exactly the same terms in g^j, multiplied by Using dehnition fl4^ . we have explicitly 
while = -2ip^ {doC)\-^-... 

This dependence of fl38|l and fl39|) on the same contraction C of gij allows us to cancel it 
in the linear combination 

_ ^2^jy00_ ( 43 ) 

On the other hand, this combination does contain a term 


- + pIb) p 


(44) 
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proportional to the second time derivative of the scalar held. Our lemmas above also prove 
that all its other terms involve at most hrst time derivatives of the helds (and up to three 
spatial derivatives, but they do not pose any difficulty for the Cauchy problem). Now, 
using the fact that the components of the Einstein tensor depend only on hrst time 
derivatives of the metric, cf. Eq. flTBD . we arrive at the following crucial combination of 
Einstein’s equations: 

(^2 ((^00 _ yOO) _ _ yaO) _ q_ ^45^ 


This generalizes Eq. (IT7|) for arbitrary n and p, and allows us to express (p in terms of 
undiherentiated helds, their spatial derivatives, and their first time derivatives. 

We now follow the same logic as in subsection III Al On taking the time derivative of 
Eq. (H5|) . we have a way to express V? in terms of helds which are diherentiated at most 
twice with respect to time. Since was the only higher-order time derivative entering 
Einstein’s equations = T"^, we conclude that all of them now become of second order 
(as far as time is concerned). Moreover, since the linear combination fl45p does not depend 
on the values n and p specifying the model C(^n+i,p), we can repeat the same argument for 
any linear combination '^np^(n,p)'^in,p)y cind this achieves our proof: Only the helds and 
their hrst time derivatives should need to be specihed on an initial value surface. 

Of course, our argument fails if the coefficient fl44|) happens to vanish at some spacetime 
point, since can no longer be expressed in terms of at most hrst time derivatives of 
helds, in such a case. This is notably what happens on a surface where pi = 0, which 


corresponds precisely to the unitary gauge chosen in Refs. |2lH23[|. However, no third time 


derivative enters any held equation when pi = 0, therefore one can keep all of them without 
any difficulty when this happens. Let us anyway mention that our argument is generic, as 
required when discussing the well-posedness of the Cauchy problem. The spacetime domains 
where Eq. (1441) may vanish have measure zero, and our proof is thus valid almost everywhere, 
in the mathematical sense of measure theory. 

Let us now show that the scalar held equation £ = 0 may also be recast as a second-order 
diherential equation with respect to time. We already saw in Eq. ([6]) that £ can be obtained 
by taking the divergence of Einstein’s equations. However, we now show how the third time 
derivatives of the metric tensor can be removed from £, while still keeping a non-zero term 
in ifi. For this purpose, consider another linear combination of Einstein’s equations 


A(G°° - T°°) + B (p„(G“° - T“°) = 0, 


(46) 


which generalizes Eq. fl?I|) above for arbitrary n and p. All second time derivatives of the 
scalar held cancel in this combination, but it still contains a term 

/ + P„B) C, (47) 

where C involves a specihc contraction of Qij, cf. Eq. (14^ . Therefore, the linear combination 
fl46p allows us to express this contraction in terms of helds diherentiated at most once with 
respect to time.® Or, by taking the time derivative of (14^ . one may express the contraction 


® Note that the same coefficient (y°A -I- enters both Eqs. (1441) and (l47l) . At the generic spacetime 

points where it does not vanish, the two linear combinations of Einstein’s equations (|45l) and (I46p can 
thus both be used, to extract ip and the specific contraction C of gij in terms of fields differentiated at 
most once with respect to time. 
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(H2|) . with Qij replaced by 9 ij, in terms of at most second time derivatives. Returning to 
the scalar field equation = 0, it turns out that the third time derivatives entering it take 
precisely the same form, namely: 


i'srd t der. — <?(<? (2n,4p+3,4p+4)9^Atl^/i4p+5 M4p+6 


'^(2n,4p-l,4p)9^Pl9^P2'^(p-l)‘^((?) ^ P4p+l P4p+2 


(48) 


As discussed below Eq. 0421) . this is a consequence of the diffeomorphism invariance of the 
theory (and of the absence of third time derivatives of the metric in T"^). Therefore, all 
third time derivatives entering £ = 0 are exactly canceled by adding to it the time derivative 
of Eq. fl46l) multiplied by 2/{(p^A + ^PaB). 

This procedure can also be extended to an arbitrary sum of Lagrangians J2np^{n,p)'^{n,p), 
although this is less obvious than for the Einstein equations, cf. the paragraph below Eq. fITSD . 
Indeed, since the coefficients A and B of Eqs. flTUD and fITT]) depend on n and p, the linear 
combination fITBD is thus specific to a single case. However, if we denote as i?(n,p) and 

C(n,p) the coefficients (IT0D -flT2 |) corresponding to a given Lagrangian C{n,p), the above results 
immediately show that the linear combination 


k / \ n,p 


kin,p)Bin,p) j V^a(G“° - T“°) = 0 
does not contain any (f, whereas its terms involving gij are of the form 




^ I ^ ^ k[n,p)^{n,p) I + I 'y ^ ^(n,p)-B(n,p) 


n,p 


n,p 


Ti^p 


^{n,p)C(^fi,p) 


(49) 


(60) 


Moreover, all its other terms involve at most first time derivatives. On the other hand, 
the third time derivatives entering the scalar field equation £ = 0 are exactly the same as 
those of —2.p>^Ylinpk{n,p)doC(n,p)) again for the same reasons as explained below Eq. (IT2D . 
Therefore, it suffices to add to this scalar field equation the time derivative of Eq. flT9l) . 
multiplied by 2 and divided by the large coefficient within the square brackets of fl50l) . to 
get a second-order differential equation with respect to time. 

Together with our previous proof that the Einstein equations themselves can be recast in 
terms of at most second time derivatives, we thus arrive at the following powerful result: 

All generalized Galileon models J2npk{n,p)£'{n,p) in curved spacetime, with arbitrary 

constant coefficients k(^ri,p), yield -|- 1 field equations which can be combined so 

that none of them involve more than second time derivatives. 


Obviously, this does not prove that all these theories are stable. For instance, it suffices 
to choose the wrong signs for some of the coefficients fc(n,p) to get a ghost scalar degree of 
freedom, notably in the simplest case +C( 2 p) with a positive sign. But this shows that the 
specific combination X)p=o'' ^(n,p)B(n,p) of Ref. [8| is not safer, nor worse, than any other 
linear combination. The only difference is that the field equations of these “covariantized” 
Galileons ([T]) do not involve any third derivative, even purely spatial, whereas we proved 
that arbitrary sums Yhnp k{n,p)B(n,p) can be cured from their third time derivatives, but they 
keep other types of third derivatives (spatial, or mixing space and time). 
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2. Horndeski/k-essence-like generalizations 


It was shown in Ref. (see also the review js^), that the above models can be generalized 
further. Indeed, the flat-space Lagrangians with p = 0, can hrst be integrated by 

parts to be rewritten as —pi ‘^( 2 ^- 2 )^^"then replace the 
factor by any function f{(p, Px) without changing the structure of the higher derivatives 
of the model. Hence the simplest case of >C( 2 ,o) = generates all k-essence theories f{ip, ip\), 
which also include the tadpole /l(i,o) = ^ and the cosmological constant /l(o,o) = const. The 
p 7 ^ 0 cases may also be generalized in the same way [^, as 


- A^^Z-2)nipA,)- ( 51 ) 

In D = 4 dimensions, one recovers the full class of Horndeski’s theories [l| , where the func¬ 
tions f(^n+i,p) of (p and ip\ multiplying the different terms need to have specihc relations 
amongst themselves in order to avoid the appearance of third derivatives in the field emo¬ 


tions. Another possible generalization of the p 7 ^ 0 cases has been considered in Refs. |21 
and is obtained by the following set of Lagrangians® 


22 


f(n+l,p}{'pi'px)^{n+l,p) f{n+l,p){,‘p-i‘px)‘A.(2n)‘ppi‘ppL2'^{p)S{q)-i 


(52) 


which are just obtained by multiplying the Lagrangians of Eq. fl26p by arbitrary functions 


/(n+i,p). The novelty claimed by Refs. |2lN23| is that even without any relation between such 


functions, these theories in any case do not generate any second scalar degree of freedom, 
although their held equations do involve third derivatives, and notably third time derivatives. 
These references concluded that the most general model in D = 4 dimensions depends on six 
arbitrary functions of <p and <p|. However, as mentioned in our Introduction, the arguments 


given in Refs. |2ll-l23l| do not appear to be fully convincing. 


Our reasoning above can be repeated without much change for these wide classes of 
theories. Indeed, let us hrst consider a single Lagrangian of the form fl5^ . It is easy to see 
that the energy-momentum tensor of such a theory is just obtained by Eq. fIM]) where one 
replaces everywhere by f(n+i,p)P>pL 2 -i and to which one adds the extra term 


- (53) 

where X = (p|, and C(n+i,p) is the generalized Galileon Lagrangian fl2^ (corresponding 
to / = 1). The three lemmas derived below Eq. fl35|) follow then immediately: No third 
derivative of the metric enters T“^, and neither T®® nor T"® contain'^ any '<P or (pi. For 
the same reasons as explained below Eq. fl42|) . one also hnds that all components of the 
stress-energy tensor involve the same contraction of Qij, namely 


T 


qO I 


'Qij 


= 


J ^ ^in+l,p) 


9ij 


(54) 


® Note that the set of theories defined respectively by Lagrangians and (l52|) can be related to each 
other by suitable integration by parts and identities between Levi-Civita tensors of the kind discussed 


^ On the other hand, there do remain some (pi in (paT°‘^ when considering a Lagrangian of the form m- 
However, even in this case, our reasoning also shows that the linear combination (1451) does not contain 
any second time derivative but those of the scalar field. 


14 












where C is given in Eq. fH2|) . Therefore, these second time derivatives of the metric still 
cancel in the same linear combination (H5|) as above, while it is easy to check that the 
coefficient of (p remains generically non zero (although it is changed with respect to the 
/ = 1 case). In conclusion, Eq. fH5]) still allows us to express (p in terms of undifferentiated 
fields, their spatial derivatives, and their first time derivatives, and the time derivative of this 
Eq. fl45p can thus be used to remove all third time derivatives entering Einstein’s equations. 
Moreover, since the linear combination fl45p does not depend on the precise model fl52p under 
consideration, one can obviously repeat the same argument for any sum of such Lagrangians. 

Let us now consider the scalar field equation S = 0 deriving from a Lagrangian ([52]). It 
involves third time derivatives of the metric, of the form 


S\;- = 




(56) 


They are proportional to those entering Eq. (IHTp . and this is not a surprise since this 

is again a consequence of diffeomorphism invariance, as explained below Eq. fl42ll . It thus 
suffices to construct the unique linear combination of the Einstein equations = 0 

and — T"°) = 0 such that all cancel, and the time derivative of this combination 

allows us to replace all third derivatives entering = 0 in terms of at most second time 
derivatives. This conclusion can also be extended to an arbitrary sum of Lagrangians (l5^ . 
by following the same reasoning as in Eqs. (Hop - flSUP above. [Let us also recall that for 
any theory, the scalar field equation S = 0 can always be recovered from the divergence of 
Einstein’s equations, Eq. (|6|).] 

In conclusion, any sum of Lagrangians J2np f{n,p) he., the k-essence-like 


extensions of generalized Galileon theories fl26l) . yield field equations which can be recast 
as a set of -p 1 second-order differential equations (as far as time is concerned). In 

D dimensions, there are + Ij classes of models, each depending on an arbitrary 


' <- Z J I- Z -> ' -*—I I 1 

function /(n,p)((p, (p|). This is consistent with the 6 classes found in Refs. j2ll-l23| for D = A, 
and would give for instance 30 classes of models in the D = 10 dimensions of string theory. 


III. HAMILTONIAN ANALYSIS OF THE QUARTIC GALILEON 


In this final section, we work m D = A dimensions and focus solely on the Lagrangian 
'^(4,0)- Our aim is to carry out a Hamiltonian analysis of this particular theory without 
fixing any gauge, and to show — from a Hamiltonian point of view — that it cannot contain 
4 degrees of freedom. As such, the results of this section support the results of the previous 
section, though the approach is different. 


A Hamiltonian analysis of “beyond Horndeski” theories was carried out in j2ll-l23|. though 


in that paper the authors restricted their attention to the unitary gauge, t = p>. As we shall 
see below, in the unitary gauge the Hamiltonian analysis is greatly simplified. Indeed, in 
an arbitrary gauge, the Lagrangian (from which the Hamiltonian is constructed) explicitly 
contains a term in (p multiplied by first time derivatives of the metric: this is the term 
which generates third time derivatives in the equations of motion. However, this term is 
also multiplied by spatial derivatives of ip which vanish in the unitary gauge. Hence in the 
unitary gauge the Lagrangian contains no second-order time derivatives. Their presence (in 
a gauge-invariant calculation, as discussed here) renders a Hamiltonian analysis much more 
involved as we shall see. 
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( 56 ) 


In this section we carry ont a Hamiltonian analysis for the non-gange fixed action 

S = J [R + /I(4,o)] , 

where we recall that /l( 4 ,o) was given in Eq. ([7]). However, precisely because of the presence 
of the ip terms mentioned above, the calculation will quickly become complex and lead to 
very sizeable expressions. Thus for reasons of clarity, we will not display all expressions in 
their full gory detail: a courageous reader is referred to Appendix [B] for more details. In 
fact, we shall push the calculation only as far as to be able to conclude that the theory 
necessarily possesses less than 4 Lagrangian degrees of freedom. 


A. ADM parametrization and primary constraints 


In the ADM parametrization (with lapse N, shift A*, and spatial metric 'jij), action (l56|) 
becomes 




where 


dtd^xN^{KijK^^ -K^ + 

dtd X ^ 6^ 6 /j; 1^ (p ‘2^i^iSQ()Sjrfi “h ^.ipiip^SQrn^Oj 

dtd^x [2(p(peSi^Sjn - ApiPeSomSjn] 

/ N NP\ 
dtd^xN^ ( 1 - 


d^xL, 




(57) 

(58) 

(59) 


and Kij is the extrinsic curvature (see Appendix lAll . 

The hrst line, linear in the lapse, is the usual ADM decomposition of the Einstein Hilbert 
action in General Relativity (GR). Notice, however, that due to the other terms, the La¬ 
grangian L is no longer linear either in N nor A*. Furthermore, L generically contains 
products of second (covariant) derivatives of p (the s^^,), which in turn contain Ghristoffel 
symbols — which are expressed in terms of time derivatives of the lapse and shifts, see 
Appendix 1^ Thus the action depends explicitly and non-linearly on A and A*, as opposed 
to the the case of GR. Notice that the term oc PiPtSQQSjm (mentioned above) in principle 
generates 3rd order derivatives in the equations of motion, though it vanishes in the unitary 
gauge. 

More generally, the non-linear dependence of action S on the variables makes it 
very difficult to invert in terms of its conjugate momenta = SL/Sjij. To alleviate 
this problem, we proceed by linearizing the action in second derivatives. [Note that this 
technique would also work for higher-order Galileon theories, for instance 71(5_o)-] That is, 
we rewrite the action fl58|) as 






d^x (s^, - V^V.v?), 


(60) 
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where and are just symmetric tensors considered as dynamical fields. The field 
is a Lagrange multiplier imposing the relation ([59]). It is straightforward to check that the 
equations of motion following from the two actions fl58|) and fl60|l are equivalent. Notice, 
however, that the price to pay for this linearization is the introduction of new degrees of 
freedom: indeed the dynamical fields are now 

( 61 ) 

where, for computational simplicity we choose to work with 

X^’' = (62) 


rather than X^". Thus there are a total of 31 dynamical fields, together with their 31 
conjugate momenta defined by 


_ 6L 

'^N = 777 , 

6N 
_ 5L 


TTi = 


= 
flU - 


SL 

6 N^ ’ 
5L 


= 


5L 


5% 




vrr = 


5L 


5s 




(63) 


The fields and their conjugate momenta satisfy the standard Poisson-Bracket (PB) relations, 
for instance {N[x),'Kis[{y)} = 6^{x,y) (see Appendix iBl for the remaining — obvious — 
relations). 

The canonical momenta are determined directly from fl57|l and fl60|l . and we find 


TT' 


pi/ _ 


(S) 


= 0 


TT, 


(A) _ 


Oi 


= 0 , 




Try ' = 0, 


TTj = 


TT^ = - (A°V - iVV,) , (64) 


where 

(65) 

The remaining conjugate momenta vr,^ and tt*-^ are given in Appendix (B] equations (1B4I1 
and (1B511 respectively. In particular depends linearly on both A°° and N [actually on 
cIo(A°°A^)]. Hence, given the expression for tt^ in fl65|) . naively one might expect this theory 
to contain two scalar degrees of freedom, A°° and ip. The momentum vr*-^ conjugate to jij 
depends on 'yij (as expected) as well as ip. 

The 23 relations in (l6T|) define 23 primary constraints 


(«) 


- ^(S) 


(A) - ^(A) 


$ 


'^Oi 


0 , 

0 , 

0 , 

00 , 


= TTi - X^'^ifi ^ 0, 

A°° 

4>7v = tttv - — (tt^ - N^ifi) ^ 0. 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


It is straightforward to see that on shell, all the primary constraints have vanishing Poisson 
brackets amongst each other. For the following discussion, it will be useful to denote the 
primary constraints by 








P= 1, 


,18. 


( 71 ) 
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That is we separate out <hj and from the remaining <hp primary constraints, since, 
as we shall see below, they must be considered differently from the others. Recall that in 
GR, and correspond, respectively, to the primary constraints vTj 0, ttat 0 which 
constitute 4 of the 8 hrst class constraints associated with diffeomorphism invariance. In a 
similar way, 4 of the primary constraints in 071 p (namely, a linear combination of and 
^isf with the other primary constraints) must also be hrst-class constraints. There must 
also be 4 secondary constraints of first class — which, in turn, do not generate tertiary or 
higher generation constraints — as a consequence of the fact that diffeomorphism invariance 
is expressed inhnitesimally through four independent paramete rs ( the four components of a 
vector) which are just differentiated once with respect to time [la, 1^ ■ 


B. Canonical Hamiltonian and secondary constraints 


The remaining Hamiltonian analysis follows the standard route (see for instance [37|]), but 
is rather involved due to the large number of fields and the intrinsically non-linear nature of 
the problem. The hrst step is to calculate the canonical Hamiltonian He from which, as a 
second step we determine the secondary constraints, obtained by imposing the preservation 
of primary constraints under time evolution. We hnd 


Hr = 


- K^) - + 2{D,Nj)H^ - + ( tt ^ tt ^) 

-(Ptt^ + g)7r^ - V 

+ + 2A°*W) - + iV(4iV)] 


—TTat 


iV"(4iV) + 


(72) 


where P, q and V are various combinations of functions appearing in the action (]57p . as 
dehned below. To simplify expressions, we dehne 


^iijm _ 


e k 


so that 


P = 

q = 
V = 


iijr 


N 






N 


and the extrinsic curvature Kij is given in terms of the canonical momenta by 




where 


TTn 




2N 




A»» 


A»» 


- (A“oy> + A“iviy‘>) 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 
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Notice that, because of the term linear in the canonical Hamiltonian appears at hrst 
sight not to be bounded from below, a la Ostrogradski [T^. However, we shall see that vr^ 
is in fact not independent of the other helds, because of one of the secondary constraints (to 
be precise it is the constraint "Hq) dehned in flSB]) and given explicitly in equation (IB8D in 
Appendix[B]). As a result the canonical Hamiltonian will actually vanish on shell, as expected 
for a diffeomorphism-invariant theory, see Eq. fl87p . Following the standard procedure, the 
total Hamiltonian is then given by 


Ht = 



(s) 


+ 8(A) 


+ C(%<1>17 + <% + (wt-N 


(79) 


thus introducing 23 Lagrange multipliers, (. 

Imposing the conservation of the primary constraints, schematically $(x) = {^{x), Hx}, 
enables us to determine 23 corresponding secondary constraints (denoted by k = <I)((r)). 
Using Eq. fl72|) together with the primary constraints in fITOl) . we hnd 


Kg) = (80) 

Kg) = 4^ + 2A°*, (81) 

Kg) = AP'' - 

{2N,so„, + {N^ - NfN^) s,m} , (82) 


(A) 

4 


K<)> = 2s„, - 2(D,Di,p)N^ - 2Nd, j + 




DiDjip 




Aoo. 


2K, 






The last 4 secondary constraints, denoted by 


(83) 

(84) 


H„ = -{4>„,i/T}, W. = -{$., ffr}. (85) 

in analogy with the Hamiltonian and momentum constraints in GR, are given in Appendix 
[B1 The hrst, "Ho contains a term in while the second 'Hi contains a term in DjTi\. In 

linear combinations with the other (primary and secondary) constraints, they must therefore 
constitute the 4 remaining secondary hrst-class generators associated with diffeomorphism 
invariance. 

In analogy with the primary constraints, it will be useful to write the set of secondary 
constraints as 

P = l,...,18. (86) 

Finally, in terms of the primary and secondary constraints, the canonical Hamiltonian in 
Eq. fl72]) can be expressed as 


H, = I d^x\NHo + N{W - (2Kg)Soo + - 4^A*^' 


4> 


N 


+ 2 ^^ + — j + ^a.(A'A“) + N%N 


(87) 


That is, as expected for a diffeomorphism-invariant theory, the Hamiltonian vanishes on 

shell. 
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C. Counting degrees of freedom, first and second-class constraints 

The theory has 2 x 31 = 62 Hamiltonian degrees of freedom, and so far we have identified 
23 primary and 23 secondary constraints. Of these, at least 8 must be hrst class (due to 
diffeomorphism invariance). If all the remaining constraints were second class and z/there 
were no tertiary constraints, then at this stage we would conclude that the theory contains 

62 - (2 X 8 ) - (46 - 8 ) = 8 ( 88 ) 

Hamiltonian degrees of freedom. That is 4 Lagrangian degrees of freedom: 2 for the graviton 
and 2 scalars. 

The next step therefore consists in determining whether or not the remaining primary 
and secondary constraints are of second class. To do so, as per the standard procedure, one 
must calculate the 46 x 46 matrix of their Poisson brackets. We order the constraints with 
the primary constraints first followed by the secondary ones. At this stage it is useful to 
notice that this 46 x 46 matrix is of the anti-diagonal form 

J e)- 

where first 23 x 23 block vanishes since the primary constraints all commute (see above), 
and the 23 x 23 block A is itself symmetric because of the Jacobi identity: 

n,} = Ht}} = {Ht, 4>.}} - {Ht, {4>„ 4>,}} = +{4>„ (90) 

We now focus on a subset of the primary and secondary constraints, namely the $p 
and Kp] the remaining constraints will be discussed later. The Poisson brackets of these 
constraints dehne a 36 x 36 matrix of the same anti-diagonal form, where the antidiagonal 
block is the 18x18 matrix 

A4i8 = {$p,«:q}, (P,Q) = 1,...,18. (91) 

The determinant of this 36 x 36 matrix is = —(det(A4i8))^, and using the commutators given 
in appendix [B] one can check that det(AJi 8 ) 7 ^ 0. Hence these 18 x 2 = 36 constraints are 
necessarily second class. This also implies that they do not generate any tertiary constraints, 
because {k, Ht} ~ 0 actually fix the corresponding Lagrange multipliers ( entering the total 
Hamiltonian in Eq. 0791) . 

However, as soon as we consider the 19x 19 matrix consisting of Mis but augmented by the 
further primary constraint 4)°°^, and the further secondary constraint nlfd, the determinant 
vanishes. Hence there exists a linear combination of and $p which commutes with 
all the ($°°^, <l)p, Kp). The explicit expression for this linear combination, denoted by 
k“). is given in Appendix [B], equation 0B16p . Similarly we have shown that there exists a 
linear combination, denoted by Kqq, of Kqq with the Kp’s and the primary constraints which 
commutes with all the (‘h(’°), d’p, , Kp). 

At this stage, the Hamiltonian analysis of £4 must necessarily proceed in one of the two 
following ways: 
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1. Either and also commute with the remaining constraints (namely "Ho, 

and <hi), and hence constitute two more first-class constraints in addition to the 
2x4 implied by diffeomorphism invariance. In this case the theory has 6 Hamiltonian 
degrees of freedom (2 for the graviton and 1 for the scalar held). 

2. Or and Rqq do not commute with the remaining constraints (namely "Hq) 'hti, 

<hjv and <hj). In that case there exist a tertiary and possibly quaternary (or higher) 
constraints. Indeed, the vanishing of the above 19 x 19 determinant implies that the 
Lagrange multiplier entering fl79|l is not hxed by the conservation of the secondary 
constraint, ~ 0. 

Whichever is correct, we arrive at the important conclusion that there are necessarily less 
than 8 Hamiltonian degrees of freedom. 

In order to see whether the hrst option is the true one, we note that if they are indeed 
hrst-class generators, and must generate a new symmetry. However, we have de¬ 
termined the transformations of the helds induced by these generators, and found that the 
action of the theory in not invariant. Hence, necessarily we arrive at the conclusion that 
case 2 is correct: there must exist a tertiary (and even perhaps a quaternary) constraint 
beyond the primary and secondary constraints found above. These are generated by the 
conservation of Kqq . Unfortunately, though, because of the size of the corresponding Pois¬ 
son brackets, the computation of these tertiary/quaternary constraints is a cumbersome task 
we felt unnecessary to attempt. 

However, we have carried out an identical calculation in the simpler case of >C( 2 ,o) (using 
the same technique with our Lagrange multipliers fl60|l . although of course this greatly 
complicates the calculation in this very simple case). Here we have shown that there are 
indeed a tertiary constraint generated by the conservation of Kqq, and even a quaternary 
constraint generated by the conservation of this tertiary one. There are no higher order 
constraints because the conservation of the quaternary constraints actually imposes the 
value of the Lagrange multiplier Qq . It happens that the secondary and tertiary do not 
commute with each other, and that the primary and the quaternary also do not commute 
with each other. Therefore, all of these four constraints are actually of second class. As 
compared to Eq. flHHj) . there are thus 2 extra constraints, which give a hnal number of 
6 Hamiltonian degrees of freedom. That is, 3 Lagrangian degrees of freedom: 2 for the 
graviton and a single scalar. We expect the same result to hold true also for /l( 4 ,o), although 
we did not do this explicit calculation because of its complexity, and due to the conclusions 
of Sec. HI] which we felt rendered it unnecessary. 


IV. CONCLUSIONS 

In this work, we have considered Galileon models in D dimensions, as well as the models 
dehned by the counterterms which have been introduced to maintain the second-order nature 
of the held equations when both the metric and the scalar are made dynamical. We have hrst 
shown that in one given such model, all the third time derivatives which appear in the held 
equations can be eliminated, leaving a set of -|- 1 held equations with at most second 

time derivatives of the dynamical helds. The same has been shown to hold for an arbitrary 
linear combination of such models, as well as their k-essence-like generalizations involving 
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free functions of cp and cpI- (In D dimensions, these models can depend on ["f + IJ 


independent such functions.) This supports the claim made previously [2l|, l22| that the 
number of degrees of freedom in these theories is only 3, counting 2 for the graviton and 
1 for the scalar. However, it does not provide an absolutely rigorous proof of this claim 
which would require a detailed Hamiltonian analysis of these models. Such a Hamiltonian 
analysis has been carried out in the second part of this paper for one of the models under 
consideration, reaching the conclusion that the number of degrees of freedom is indeed 
strictly less than 4. It seems, however, very likely, in light of our results, and taking into 
account the gauge invariance of the theory (which has been fully kept in our analysis, in 
contrast to previous works) that the final number of degrees of freedom is well only 3. Indeed, 
using the first part of this paper, one expects that the “reduced” field equations (containing 
only second derivatives) can still be decomposed into 4 Lagrangian constraints (coming 
from the invariance of the theory under reparametrization) plus 6 dynamical equations 
for 6 dynamical metric variable and an other one for the scalar. Note, however, that the 
extraction of these would be 4 Lagrangian constraints is much trickier than in standard 
general relativity because it can be checked that, in an ADM language, the field equations 
contain second time derivatives of the spatial part of the metric, but also of the lapse and 
the shift. Provided these Lagrangian constraints exist, the gauge invariance would then 
reduce to two the number of dynamical components in the metric. Rigorously speaking, our 
Hamiltonian analysis also leaves open the possibility that there is just one tertiary second 
class constraint (and no quaternary). This would result in an odd (times inhnity®) number 
of second class constraints, which can happen in a held theory (see e.g. jssf and references 
therein), but seems to us unlikely in a bosonic and Lorentz-invariant theory. This, however, 
deserves further investigation to be rigorously checked. 
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Appendix A: Christoffel symbols in ADM variables 


In the ADM notation, the 4D metric and its inverse are written as 


f-N^ + NkN^ Ni\ 
V Tj / 


9 ^'^ = 


1 \ 
Ar2 \ 

ij _ j ’ 

7 7V2 / 


(Al) 


In classical mechanics, there is always an even number of second-class constraints, but this is no longer 
true in continuous field theories, where each constraint should actually be understood as an infinity of 
them, since it is imposed at every point of the Cauchy surface. 
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where here and in the following all spatial indices i,j, k,... are raised and lowered with the 
spatial metric 7 *^. On nsing (lAlj) the Christoffel symbols can be calcnlated. These will 
contain 7 ^, N and iVh The covariant derivative associated with the metric 7 *^ is denoted 
by D, and Kij is dehned by 

^ij ~ (77 ~ DiNj — DjNi). (A 2 ) 

We hnd 

Too = ^ (^ + , 

rSo = +1^‘ + 2NN<K,t h’’-' - + iVV‘(B,iVO 

„ / ,1- N>N‘‘\ 

+Nd,N , 

r»„ = 1 (diN + K„N ^), 

N^d-N 1 

T'o =- + + 


■pO 

'-ji 




-pn 


■^^T + r^z(7), 


(A3) 


where r”;( 7 ) stands for the Christoffel symbols of the spatial metric 'jij. 


Appendix B: Hamiltonian analysis of ^^(4,0) 


We give here some of the intermediate expressions reqnired for Sec. IIIII We will use one 
abuse of notation in this Appendix (and only here): namely we denote Y^(pj by (p\ 

The 31 dynamical variables given in fl6T|) . namely 


iV,iV\7p 

5 ■) 

(Bl) 

and their conjugate momenta 



'^N 1 '^ 5 


(B2) 

satisfy the Poisson Brackets 



{A(x),7r7v(l/)} = S^{x,y), 
{N\x),7rj{y)} = 6}6^{x,y), 
Hj{x),7rP‘^{y)} = (5j)hj)h^(x,i/), 
{A°°(x),7r^o(l/)} = 

{soo(x),7r^°(?/)} = S^{x,y), 


{X^\x),Tioj{y)} = 5j5^{x,y), 

{So*(a;), 7^5^( 1 /)} = Sj6^{x,y) , 


{A*^(a:),<(i/)} = 6i:6{^6^{x,y) , 

{sij{x),T^l\y)] = hj’.5j)5^(x,|/) , 

(B3) 
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with all other commutators vanishing. 

Starting from fl57|) . fl58|) and fl60|) . the conjugate momenta 'K^p and are found to be 
given by 


Tiw — + Q + + 


A°0iV , KijQ 


N 


+ 


N 


= ^(K» - A'7«) + + 1^. 


where 


Q = ^(5.]V)VA“> + (AfA“) , 

Qij = + 2A°W^'+ A®^ 


rpimn _ 2iViV’®A°° ( 7®"® — 

On using 


2iV2 


+ 2A°®®A^ ( 7®®” - 


iVW® 

Ar2 


-A® 


^iV® 


, d70|) and fl72|) . one can calculate the secondary constraints 

•Ho = 

H, = 


The first is given by 

Hn 


^ [(A-yA'-J - le) - P>A] + K„B‘= + ^ + ^ 


^ [P(ir")‘" + (-n + U) - (AVi) (2!r'P + <))] 


+2A“a. + a, (]VA“V‘ + irjvV) + ^N‘{dtN) 


N 


Aoo 


+ —^®iV^(D,iV0 - (a.TT^) 


/A°°iV® 




iV 


$ 


TV 


N 


K, 




iVTT. 


A^h/VT pi \ 1 

iV®iV^ + 2^^ + — +iV^4iV + 2—a.(iVA°®)+ 


where 


N 


Bij = ^ ( 2 N^N^ + 2 


AO(ijYi) yj 


+ ^) + 2A“>Af‘VI, 


(B4) 

(B5) 


(B6) 


(B7) 


,(B8) 


(B9) 


and P and q were defined in (17^ and (jTHjl respectively. The quantities Vi and V 2 are 
different components (depending on their A^-dependence) of V defined in Eq. fl76|) b Namely 
V = Vi + V 2 with 


VA = {somSoj - SooSjm) - N^So^S 

V 2 = 
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(Bll) 



















Finally, 


Hi{x) = -2DjTT\-2NK, 






N 


XOj 


Aoo 


+ 2Dj (A°°ArVi) + 2{DjD,ip)X^^ - - ttnO^N + 

+ (PiTT^ - Pi{2Pn^ + q) 

I Jmn, 


N 


[Ni {e’^^'^SjraSknPq) “ e"^•S.■n (TT-'S 


''Fq. 

2'^jn ^qm 


2^qSom,)] • 


(B12) 


We now list the different commutators of primary and secondary constraints, which can 
be straightforwardly determined from the primary constraints given in (l66ll - (17ni) and the 
secondary constraints given in fl5UD - fl5^ . In particular, we hnd that satishes 




N 


(B13) 


with all its remaining commutators (including with the vanishing. Then we hnd 


= - 26 ], 

{■i’i'A't™} = 0, 

{*S). '*” } = ^ - 6«<!'6»>'V/43((iV2 - iV,]V»)] , 

(B14) 


where this last expression should be understood to be symmetrized over both pq and ij. 
Finally 


{Ofib'tg’} 

r<T,F) 

l^Oi ; ’^mnS 

5 ^mnJ 


N 


[{^) 


VTjv 


2^Aoo 


TTat 


(iVAOO 




[X'mi'ynf T 'yni'Jfm ; 
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N 


AN^iX' 


00'|2 


[XmiXnj T 'Xmj^ni 'yijXmn] 


(B15) 


From the 19 x 19 matrix discussed in Sec. IIII Cl the following linear combination of 
primary constraints (denoted by commutes with all the {d’p, np}: 


1.00 = [N^ - 2iVV?vr^(vr^ - N’^Pk) + vr^(vr^ - N’^Pk){Tr^ - 2N%)] <Fg 

- [N^pyiTT^ - ivVfc) + N^Pj{Nyy{7r^ - Nyf] pyi^ 

- (P - JVVOVi-l’Sii'.’i - 2iV(¥.‘<*!"V - (B16) 
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We have also proved that it commutes with <hAr as well as "Ho- The secondary constraint Kqq 
which commutes with all the {*hp} as well as <hAr is given by the same expression but where, 
on the right hand side, the primary constraints are replaced by secondary constraints k 
(with the same labels). One must also complement this with a linear combination of 
primary constraints so that it commntes with all secondary constraints [see onr discussion 
below Eq. CT ]. 
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